Linearized Navier-Stokes Problem Research Articles

An Eulerian finite element method for the linearized Navier–Stokes problem in an evolving domain

Abstract The paper addresses an error analysis of an Eulerian finite element method used for solving a linearized Navier–Stokes problem in a time-dependent domain. In this study, the domain’s evolution is assumed to be known and independent of the solution to the problem at hand. The numerical method employed in the study combines a standard backward differentiation formula-type time-stepping procedure with a geometrically unfitted finite element discretization technique. Additionally, Nitsche’s method is utilized to enforce the boundary conditions. The paper presents a convergence estimate for several velocity–pressure elements that are inf-sup stable. The estimate demonstrates optimal order convergence in the energy norm for the velocity component and a scaled $L^{2}(H^{1})$-type norm for the pressure component.

A Two-Level Nonconforming Rotated Quadrilateral Finite Element Method for the Stationary Navier–Stokes Equations

In this paper, we propose a two-level nonconforming rotated finite element (TNRFE) method for solving the Navier–Stokes equations. A new nonconforming rotated finite element (NRFE) method was proposed by Douglas added by conforming bubbles to velocity and discontinuous piecewise constant to the pressure on quadrilateral elements possessing favorable stability properties. The TNRFE method involves solving a small Navier–Stokes problem on a coarse mesh with mesh size H and a large linearized Navier–Stokes problem on a fine mesh with mesh size h by the NRFE method. If we choose h = O H 2 , the TNRFE method gives the convergence rate of the same order as that of the NRFE method. Compared with the NRFE method, the TNRFE method can save a large amount of CPU time. In this paper, the stability of the approximate solutions and the error estimates are proved. Finally, the numerical experiments are given, and results indicate that the method is practicable and effective.

Augmented Lagrangian preconditioner for large-scale hydrodynamic stability analysis

Hydrodynamic linear stability analysis of large-scale three-dimensional configurations is usually performed with a“time-stepping” approach, based on the adaptation of existing solvers for the unsteady incompressible Navier–Stokes equations. We propose instead to solve the nonlinear steady equations with the Newton method and to determine the largest growth-rate eigenmodes of the linearized equations using a shift-and-invert spectral transformation and a Krylov–Schur algorithm. The solution of the shifted linearized Navier–Stokes problem, which is the bottleneck of this approach, is computed via an iterative Krylov subspace solver preconditioned by the modified augmented Lagrangian (mAL) preconditioner (Benzi et al., 2011). The well-known efficiency of this preconditioned iterative strategy for solving the real linearized steady-state equations is assessed here for the complex shifted linearized equations. The effect of various numerical and physical parameters is investigated numerically on a two-dimensional flow configuration, confirming the reduced number of iterations over state-of-the-art steady-state and time-stepping-based preconditioners. A parallel implementation of the steady Navier–Stokes and eigenvalue solvers, developed in the FreeFEM language, suitably interfaced with the PETSc/SLEPc libraries, is described and made openly available to tackle three-dimensional flow configurations. Its application on a small-scale three-dimensional problem shows the good performance of this iterative approach over a direct LU factorization strategy, in regards of memory and computational time. On a large-scale three-dimensional problem with 75 million unknowns, a 80% parallel efficiency on 256 up to 2048 processes is obtained.

Parallel two-grid finite element method for the time-dependent natural convection problem with non-smooth initial data

In this paper, we consider the parallel two-grid finite element method for the transient natural convection problem with non-smooth initial data. Our numerical scheme involves solving a nonlinear natural convection problem on the coarse grid and solving a linear natural convection problem on the fine grid. The linear natural convection problem can be split into two subproblems which can be solved in parallel: a linearized Navier–Stokes problem and a linear parabolic problem. We firstly provide the stability and convergence of standard Galerkin finite element method with non-smooth initial data. Secondly, we develop optimal error estimates of two-grid finite element method for velocity and temperature in H1-norm and for pressure in L2-norm. In order to overcome the difficulty posed by the loss of regularity, some suitable weight functions are introduced in our stability and convergence analysis for the natural convection equations. Finally, some numerical results are presented to verify the established theoretical results.

Two-Level Newton Iteration Methods for Navier-Stokes Type Variational Inequality Problem

AbstractThis paper deals with the two-level Newton iteration method based on the pressure projection stabilized finite element approximation to solve the numerical solution of the Navier-Stokes type variational inequality problem. We solve a small Navier-Stokes problem on the coarse mesh with mesh size H and solve a large linearized Navier-Stokes problem on the fine mesh with mesh size h. The error estimates derived show that if we choose h = O (|logh|1/2H3), then the two-level method we provide has the same H1 and L2 convergence orders of the velocity and the pressure as the one-level stabilized method. However, the L2 convergence order of the velocity is not consistent with that of one-level stabilized method. Finally, we give the numerical results to support the theoretical analysis.

Preconditioning the Incompressible Navier-Stokes Equations with Variable Viscosity

The aim of the project is to construct, analyse and implement fast and reliable numerical solution methods to simulate multi-phase flow, modeled by a coupled system consisting of the time-dependent Cahn-Hilliard and incompressible Navier-Stokes equations with variable viscosity and variable density. This thesis mainly discusses the efficient solution methods for the latter equations aiming at constructing preconditioners, which are numerically and computationally efficient, and robust with respect to various problem, discretization and method parameters.In this work we start by considering the stationary Navier-Stokes problem with constant viscosity. The system matrix arising from the finite element discretization of the linearized Navier-Stokes problem is nonsymmetric of saddle point form, and solving systems with it is the inner kernel of the simulations of numerous physical processes, modeled by the Navier-Stokes equations. Aiming at reducing the simulation time, in this thesis we consider iterative solution methods with efficient preconditioners. When discretized with the finite element method, both the Cahn-Hilliard equations and the stationary Navier-Stokes equations with constant viscosity give raise to linear algebraic systems with nonsymmetric matrices of two-by-two block form. In Paper I we study both problems and apply a common general framework to construct a preconditioner, based on the matrix structure. As a part of the general framework, we use the so-called element-by-element Schur complement approximation. The implementation of this approximation is rather cheap. However, the numerical experiments, provided in the paper, show that the preconditioner is not fully robust with respect to the problem and discretization parameters, in this case the viscosity and the mesh size. On the other hand, for not very convection-dominated flows, i.e., when the viscosity is not very small, this approximation does not depend on the mesh size and works efficiently. Considering the stationary Navier-Stokes equations with constant viscosity, aiming at finding a preconditioner which is fully robust to the problem and discretization parameters, in Paper II we turn to the so-called augmented Lagrangian (AL) approach, where the linear system is transformed into an equivalent one and then the transformed system is iteratively solved with the AL type preconditioner. The analysis in Paper II focuses on two issues, (1) the influence of a scalar method parameter (a stabilization constant in the AL method) on the convergence rate of the preconditioned method and (2) the choice of a matrix parameter for the AL method, which involves an approximation of the inverse of the finite element mass matrix. In Paper III we consider the stationary Navier-Stokes problem with variable viscosity. We show that the known efficient preconditioning techniques in particular, those for the AL method, derived for constant viscosity, can be straightforwardly applicable also in this case.One often used technique to solve the incompressible Navier-Stokes problem with variable density is via operator splitting, i.e., decoupling of the solutions for density, velocity and pressure. The operator splitting technique introduces an additional error, namely the splitting error, which should be also considered, together with discretization errors in space and time. Insuring the accuracy of the splitting scheme usually induces additional constrains on the size of the time-step. Aiming at fast numerical simulations and using large time-steps may require to use higher order time-discretization methods. The latter issue and its impact on the preconditioned iterative solution methods for the arising linear systems are envisioned as possible directions for future research.When modeling multi-phase flows, the Navier-Stokes equations should be considered in their full complexity, namely, the time-dependence, variable viscosity and variable density formulation. Up to the knowledge of the author, there are not many studies considering all aspects simultaneously. Issues on this topic, in particular on the construction of efficient preconditioners of the arising matrices need to be further studied.

A Comparative Study of Block Preconditioners for Incompressible Flow Problems

Problem statement: We consider the numerical solvers for the linearized Navier-Stokes problem. Both the Stokes problem and Oseen problems are considered. Approach: We used the Mark and Cell (MAC) discretization method to discretize the Navier-Stokes equations. We used preconditioned Krylov subspace methods to solve the resulting linear systems. Results: Numerical experimental results are performed to compare the different preconditioners. Conclusion: The choice of the preconditioner is highly problem dependent and we give the suggestions for individual cases.

Field-of-Values Convergence Analysis of Augmented Lagrangian Preconditioners for the Linearized Navier–Stokes Problem

We study a block triangular preconditioner for finite element approximations of the linearized Navier-Stokes equations. The preconditioner is based on the augmented Lagrangian formulation of the problem and was introduced by the authors in [SIAM J. Sci. Comput., 28 (2006), pp. 2095-2113]. In this paper we prove field-of-values type estimates for the preconditioned system which lead to optimal convergence bounds for the GMRES algorithm applied to solve the system. Two variants of the preconditioner are considered: an ideal one based on exact solves for the velocity submatrix, and a more practical variant based on block triangular approximations of the velocity submatrix.

Optimal Boundary Control of the Incompressible Fluid Problems in Rotation form

Problem Statement: We consider the optimal boundary control of the linearized Navier-Stokes problem. Both the Stokes problem and Oseen problem in rotation form are considered. Approach: We use the Mark and Cell (MAC) discretization method to discretize the optimization problem with linear constraints including the Stokes problem and the Oseen problem in rotation from. Then Reduced Hessian methods are to solve the problem. Results: Numerical experimental results are performed for the boundary optimization problem with the Stokes constraints and Oseen constraints. All the computed solutions and the desired solutions are compared. Conclusions: The proposed reduced Hessian methods have a high accuracy obtaining the optimal boundary conditioning for the Stokes problem and the Oseen problem in rotation form.

Multi-level spectral Galerkin method for the Navier–Stokes equations, II: time discretization

A fully discrete multi-level spectral Galerkin method in space–time for the two-dimensional nonstationary Navier–Stokes problem is considered. The method is a multi-scale method in which the fully nonlinear Navier–Stokes problem is only solved on the lowest-dimensional space $H_{m_{1}}$ with the largest time step Δt 1; subsequent approximations are generated on a succession of higher-dimensional spaces $H_{m_{j}}$ with small time step Δt j by solving a linearized Navier–Stokes problem about the solution on the previous level. Some error estimates are also presented for the J-level spectral Galerkin method. The scaling relations of the dimensional numbers and time mesh widths that lead to optimal accuracy of the approximate solution in H 1-norm and L 2-norm are investigated, i.e., m j∼m j−1 3/2 , Δt j∼Δt j−1 3/2 , j=2,. . .,J. We demonstrate theoretically that a fully discrete J-level spectral Galerkin method is significantly more efficient than the standard one-level spectral Galerkin method.

Multi-level spectral galerkin method for the navier-stokes problem I : spatial discretization

A multi-level spectral Galerkin method for the two-dimensional non-stationary Navier-Stokes equations is presented. The method proposed here is a multiscale method in which the fully nonlinear Navier-Stokes equations are solved only on a low-dimensional space ** subsequent approximations are generated on a succession of higher-dimensional spaces ** j=2, . . . ,J, by solving a linearized Navier-Stokes problem around the solution on the previous level. Error estimates depending on the kinematic viscosity 0<ν<1 are also presented for the J-level spectral Galerkin method. The optimal accuracy is achieved when ** We demonstrate theoretically that the J-level spectral Galerkin method is much more efficient than the standard one-level spectral Galerkin method on the highest-dimensional space **.

On the nonstationary linearized Navier–Stokes problem in domains with cylindrical outlets to infinity

The nonstationary linearized Navier–Stokes system is studied in the domain with cylindrical outlets to infinity in weighted function spaces with an exponential weight function. It is proved that under natural compatibility conditions there exists a unique solution with prescribed fluxes over the sections of outlets to infinity that exponentially tends in each outlet to the corresponding nonstationary Poiseuille flow.

Singular solutions of the stationary linearized Navier-Stokes problem for micropolar viscous liquids

A system of singular solutions of the stationary linearized Navier-Stokes problem for viscous liquids characterized by an asymmetric stress tensor is constructed.